Collective dipole oscillations in hydrogenated 

metals 



J.S.Brown 

Clarendon Laboratory, University of Oxford, Parks Rd. 0X1 3PU, UK 
Abstract 

It is observed that low-lying transitions of an interstitial hydrogen adatom on a 
metallic lattice correspond classically to dipoles oscillating at frequencies where 
band electrons typically have a low electromagnetic absorption. Such transitions 
hence give rise to an essentially unscreened time-dependent perturbation of the 
lattice potential experienced by neighbouring hydrogen states. It is shown that 
above a certain hydrogen concentration, coherently phased oscillations at adjacent 
sites may result in a dipole-dipole interaction larger than the excitation energy. This 
would give rise to a new ground state for the ensemble of coupled oscillators and a 
reduction in lattice symmetry below some phase transition temperature. An explicit 
calculation is carried out for PdH and PdD, and indirect evidence for the existence 
of the postulated phase is presented. 
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1 Introduction 



It is generally accepted that certain properties of substoichiometric transi- 
tion metal hydrides are not fully explained by existing models. In addition to 
the long-standing reverse isotope dependency of some superconducting tran- 
sition temperatures, it is not clear whether the high mobility of protons in 
certain metals is fully accounted for by phonon processes. Furthermore, neu- 
tron scattering studies [1] appear to suggest a quantum phase correlation 
between hydrogen adatoms at the femtosecond timescale, which is collabo- 
rated by EELS studies [2] indicating hydrogen derealization on Ni(lll) sur- 
faces. Finally, and perhaps most intriguingly, beam-foil experiments [3,4] have 
revealed anomalously large fusion cross-sections when low energy deuterons 
traverse foils composed of an wide range of metallic elements [5,6]. A recent 
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attempt [7] to explain the latter effect in terms of phonon exchange processes 
failed to account for the full magnitude of the observed anomaly. We were 
accordingly led to look for a further, hitherto neglected, generic effect whose 
inclusion in the model Hamiltonian might yield a more complete description 
of these systems. 



2 Hydrogen adatoms on a periodic lattice 

The many-body Hamiltonian governing the interactions of interstitial hydro- 
gen adatoms amongst themselves and with a periodic lattice of arbitrary sym- 
metry can be expressed as 

+ Vij A C R+A j,C Rj C R+A ,i'CR t i 

where c Ri creates a hydrogen adatom in state i at site R and e« is an eigen- 
value of the single-particle Schrodinger equation solved in the region of the 
site. This is invariably a high symmetry (e.g. octahedral, tetragonal) point of 
the primitive cell at which the total effective lattice potential is a local mini- 
mum. 

A denotes a member of the set of near neighbour displacements, r creates a 
phonon of energy u and polarization T. The elements of the phonon-coupling 
matrix M are in general a complicated function of all five indices. The ele- 
ments of V can be derived, in principle, from the photon propagator kernel 
in the dielectric medium of the band electrons. Repeated indices indicate an 
implicit summation over that variable. 

The site-diagonal A = elements of the hydrogen-phonon coupling M can 
be transformed into an effective self-trapping effect, whereby each hydrogen 
creates a static phonon displacement field of short range. The price we pay for 
the simplification resulting from the neglect of phonon processes is that the 
site energies q and wavefunctions ipi(r — R) become somewhat dependent on 
the spatial distribution of occupancy at neighbouring sites. 
The A^O elements of M are responsible for site-hopping. This can be repre- 
sented by a term of the form t c l R+A c R,i where t is on the order of 1 meV 
and the sum is taken over nearest neighbour A only. 

The site-diagonal (electrostatic repulsion) elements Vij j0 are of the order of a 
few eVs and effectively prohibit multiple occupancy of sites in all conceivable 
lattices. 

At small energies corresponding to transitions of the band electrons about 
the Fermi level 6f, the interaction potential V is reasonably well described 
by Thomas- Fermi screening. With the values of N(e F ) typical [8] of the met- 
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als capable of absorbing large amounts of hydrogen, neighbours are effectively 
screened from one another and the intersite interaction (i.e. A 7^ 0) is of the or- 
der of a few meV only. However, this holds only to the extent that neighbouring 
occupants both reside in their respective ground states (i.e. i(t) = j(t) = 0). 
The first excited states in the potential well experienced by hydrogen adatoms 
are typically [2,9,12] on the order of 0.1 eV above the ground state. The cor- 
responding transition frequency lies in a region where the electromagnetic 
absorption coefficient K of a metal is typically of the order of a few tens of 
cm" 1 and where the attenuation over a distance of a few A is insignificant. 
This means that, for i ^ j, the intersite interaction is not subject to the 
exponential screening factor. Accordingly: 

Vij(r,r') « y^e-^ e2N ^ r ~ r '^ (2) 

where we have replaced the displacement vector index A by the continuous 
functional dependence over real space. 

This crucial distinction between slow displacements of the hydrogen adatom 
and rapid transitions between relatively long-lived widely separated levels ap- 
pears to have been overlooked in previous models of the hydrogen-metal sys- 
tem. We accordingly adopt the hypothesis that, by taking proper account of 
the frequency-dependence of the interaction between neighbours in the model 
Hamiltonian, we may be able to account for at least some of the anomalous 
effects mentioned in the introduction above. 

Incorporating all phonon effects and static nearest neighbour repulsion into 
the state energy as in [9], the original Hamiltonian (1) becomes: 

tin R ,i + Un R {n R - 1) + tc R+Ai c R;i + V itj ^ iA ^oc R+Ad/ c RJ c R+ A,i'C R ,i (3) 

where we have made the replacement n = c^c, the effective on-site repulsion 
energy is represented by a constant U and el are the neighbour-dependent self- 
trapped energies characterized by a bandwidth on the order of a few meV. 
Retaining just the lowest order term in the multipole expansion of the Coulomb 
operator ~^hp\i we have the following explicit form for the dominant (dipole- 
dipole) contribution to the expectation value of simultaneous (i.e. photon me- 
diated) transitions between two pairs of states at sites with position vectors 
R and R' = R + A 

UJ'IVaM = A- 3 0>|i) • 0>V) - 3A- 5 (j\A.r\i)(j'\A.r'\i') (4) 

This is the first-order correction to the energy of a pair of neighbours due to 
their oscillatory dipole-dipole interaction. It is clear that a significant second- 
order effect of the dipole-dipole interaction will be to break the cubic symmetry 
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of the lattice such that an > a±, since any elongation along the dipole axis will 
tend to increase the dipole moment of the odd-parity states as well as reduce 
their (single-particle) energy relative to the other members of their respective 
multiplets. 

The lowest energy configurations of O-site dipoles in an fee host and T-site 
dipoles in a bec host are shown below: 



Fig. 1. Dipole alignment in cubic cells 




Solid circles denote metal cores and open circles denote hydrogen sites. 
With randomly distributed fractional occupancy A and complete dipole align- 
ment, the sum over all fee O-sites yields a dipole-dipole energy per hydro- 
gen of ~ — 4.4e 2 <i 2 aJ 3 A. The corresponding result for T-sites in a bec host is 
« -21.2e 2 d 2 a- 3 \. 

Whether, and to what degree such an alignment actually takes place depends 
upon the temperature and the relative magnitudes of the transition and dipole- 
dipole energies. We will describe the outlines of an approximate approach to 
this problem in the following section. 



3 Mean-field theory 



We consider the effective Hamiltonian for a pair of neighbouring sites A, B. In 
the interests of clarity and brevity we will limit ourselves to just two accessible 
states. These need not necessarily belong to the two lowest multiplets. The 
extension to a more realistic number of states is relatively straightforward. 

Keff = £o( n A,o + n B ,o) + ei(n A i + n B ,i) 

(5) 

+ aV(c Bfl C ] A1 CB,lCA,Q + 4,i c Ji,o c b,oC.4,i) 
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where V is the nearest-neighbour dipole-dipole interaction given by (4) and a, 
the off-diagonal order parameter of the model, is the thermodynamic average 
over the ensemble of the term in V summed over all sites. 
Each site can be either unoccupied or contain a hydrogen in state A or B. 
There are hence nine possible states for the two sites in this model. The 
chemical potential fx must be included, since we assume an external reservoir 
of hydrogen and treat the hydrogen occupancy fraction A as a state variable. 
The effective two-site energy levels are the eigenvalues of the following matrix: 



^e// — l^E n i,R = 

00000 

Oeo-yuOO 

e 1 -fi 

OOOeo-yuO 

2e -2/i 

e + ei - 2/x a*V 

0000 Oex-^0 

aV ei + e - 2/x 

2ei - 2/x 



By inspection one sees that the two-site partition function is 

Z = l + 2e l3 ^- €o) + 2e^- ei) + e 2/3 ^- €o) + e 2/3( ^ ei) + e~ pe + + +e~ pe - (7) 



(6) 



where f3 = (kT) 1 and 



e± = e + ei — 2 fx ± aV 



(8) 



A = 2 



e /3(M-6 ) + + e 2/3( M -e ) + e 2/3( M - ei ) + e -/3e+ + +e "^-l (9) 



a = A 



(10) 



where A is the symmetry-dependent factor resulting from sum over all sites. 
The system of equations (6-10) can be solved iteratively until self-consistency 
is obtained. 
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4 Application to PdH and PdD 



These systems have been the subject of a great deal of experimental and the- 
oretical work. Adiabatic Pd-H pair potentials that yield a spectrum of lower 
excited states agreeing well with empirical determinations of the optical tran- 
sition frequencies have been obtained from ab initio DFT calculations [9,10]. 
We accordingly used the pair potential VII of [9] to solve the Schrodinger 
equation for a single hydrogen adatom over a real space mesh of width of 
0.04 A. The lowest states were all found to be centred at the octahedral (O- 
site) symmetry point equidistant from six metal ions. Some higher states are 
centred on the tetragonal (T-site) symmetry points located at the midpoint 
between the O-sites and the metal cores, but these lie higher in energy than 
those we will be considering. The A lg singlet ground state has even parity 
(+- 1— 1-) with respect to all three axes. The first excited state is a T\ u triplet of 
odd parity with respect to one axis. The second excited state is a A 2g triplet 
of odd parity with respect to two axes. Some of our results are tabulated below: 



Isotope 


proton 


deuteron 


Occupancy A 


0.25 


1.0 


1.0 


1.0 


0.25 


1.0 


1.0 


1.0 


a y ,a z (A) 


3.94 


4.07 


4.07 


4.07 


3.94 


4.07 


4.07 


4.07 


a x (A) 


3.94 


4.07 


4.22 


4.37 


3.94 


4.07 


4.22 


4.37 


e +++ 


112 


82 


87 


100 


73 


52 


59 


73 


e-++ 


200 


151 


139 


140 


130 


95 


92 


98 


e +-+> e ++- 


200 


151 


156 


167 


130 


95 


103 


116 


e_ + _, e — + 


279 


211 


200 


200 


183 


135 


131 


137 


e+__ 


279 


211 


216 


226 


183 


135 


142 


155 


(+ + +\x\ - ++) 


0.15 


0.17 


0.20 


0.22 


0.13 


0.15 


0.18 


0.20 


-AAe 2 d 2 af\ 


-24 


-28 


-37 


-47 


-18 


-22 


-29 


-38 


^■exc 


44 


34.5 


38 


48 


28.5 


21.5 


28.5 


40.5 



The state energies are all in meV relative to the lattice potential at the O- 
site. For a\\ > a± = 4.07 A, a lattice strain energy contribution, deduced from 
the Pd-Pd force constant given by [11], has been added to all states. The 
penultimate row tabulates the dipole-dipole energy per adatom at maximum 
alignment and A = 1. The final row tabulates the corresponding excitation 
energy per adatom for the pair state e_ in (8) relative to the ground state in 
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the undistorted lattice. 

It can be seen from our results for A = 1 that lattice expansion along the dipole 
axis results in a rapid increase in dipole-dipole energy and a simultaneous de- 
crease in excitation energy. Moreover, the maximum dipole energy appears to 
closely match the excitation cost of the oscillations for a wide range of lattice 
strain values. On the basis of our data, the coherent oscillating phase would 
appear to be (just) favoured in the deuteride at low temperatures, but in view 
of the inherent uncertainty in our data, it is not considered worth applying 
the method of Section 3 to deduce a phase transition temperature. 



5 Conclusion 

We have demonstrated that conventional models of the hydrogen-lattice sys- 
tem have grossly underestimated the effective coupling between neighbouring 
hydrogen sites. We have furthermore shown that, in one metal at least, the 
neglected effect may give rise to a new low temperature phase in which all 
hydrogen adatoms oscillate coherently. The posited effect has the potential 
to explain both the anomalous shoulder in the optical phonon spectrum and 
the quantum entanglement indicated by neutron scattering and EELS studies 
and we therefore deem that further theoretical investigation is warranted. A 
further possible consequence might be an enhanced probability of fusion of 
deuterons in the collective state. Owing to the translational symmetry of the 
many-body wavefunction across lattice cells, we speculate that the dominant 
decay channel might be radiationless. 
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